Photon splitting in atomic fields 



o 
o 

(N 
> 

o 

:z; 

o 



R.N. Lee, A.L. Maslennikov, A.I. Milstein, 
V.M. Strakhovenko, and Yu.A. Tikhonov 

G.I. Budker Institute of Nuclear Physics, 
630090 Novosibirsk, Russia 



Abstract 

Photon splitting due to vacuum polarization in the electric field of an atom is consid- 
ered. We survey different theoretical approaches to the description of this nonlinear 
QED process and several attempts of its experimental observation. We present the 
results of the lowest-order perturbation theory as well as those obtained within the 
quasiclassical approximation being exact in the external field strength. The exper- 
iment where photon splitting was really observed for the first time is discussed in 
details. The results of this experiment are compared with recent theoretical estima- 
tions. 
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1. Introduction 

Virtual creation and anniliilation of electron-positron pairs is known to induce a self-action 
of an electromagnetic field, which results in such effects as coherent photon scattering and 
photon splitting in an external field, and photon-photon scattering. Among these processes 
the latter was the first one subjected to the exhaustive theoretical investigation [1] (see 
also [2] and references therein). However, this process was never observed experimentally. 
At present, coherent photon scattering in the electric field of atoms (Delbriick scattering) 
is investigated in detail both theoretically and experimentally [3-5]. The amplitudes and 
cross sections of this process were calculated in the lowest order in Za (Born approxi- 
mation) for arbitrary photon energy cu and exactly in this parameter for a; ^ m, Z|e| is 
the nucleus charge, a — e^/An — 1/137 is the fine- structure constant, e and m arc the 
electron charge and mass, ^ = c = 1. It turned out that higher orders of perturbation 
theory with respect to Za (Coulomb corrections) play an important role and drastically 
modify the cross section of Delbriick scattering on heavy atoms at high photon energy. 
The experience gained at the investigation of Delbriick scattering was extremely useful for 
the study of photon splitting. In particular, the importance of the Coulomb corrections 
in the latter was recognized. 

In the process of photon splitting, the initial photon turns in the electric field of an atom 
into two photons sharing its energy uj^. For a long time, the amplitude of this process 
was known only in the Born approximation [6,2]. For uoi <^ m, the asymptotic form of 
the Born amplitude was derived earlier in [7,8]. First successful estimate of the total cross 
section of high-energy {ui 3> m) photon splitting was made in [9]. More accurate results 
were obtained in [10] within the same approach for the total as well as for the partly 
integrated cross sections. The detailed numerical investigation based on the results of 
[6,2] was performed in [11,12]. 

Recently, an essential progress in understanding of photon splitting phenomenon was 
achieved due to the use of the quasiclassical approach. In papers [13-15] various differ- 
ential cross sections of high-energy photon splitting have been calculated exactly in the 
parameter Za. Similar to the case of Delbriick scattering, the exact cross section turns out 
to be noticeably smaller than that obtained in the Born approximation. So, the detailed 
theoretical and experimental investigation of photon splitting provides a new sensitive 
test of QED when the effect of higher-order terms of the perturbation theory with respect 
to the external field is very important. 

The observation of photon splitting is extremely hard problem due to severe background 
conditions. The significance of various competing processes depends on photon energy. For 
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instance, at oui < m the dominant background process is double Compton scattering on 
atomic electrons. The search for photon splitting in this energy region was performed in 
two experiments [16,17]. As pointed out in [16,17], the results of both experiments do not 
agree with the theoretical estimates. Nevertheless, the authors of [16] argue that their re- 
sults strongly indicate the existence of photon splitting. The energy region a;i ^ m is more 
favorable for the observation of the phenomenon, though the background is still rather 
severe. Using the advantage of the intense source of tagged photons, the first successful 
observation of photon splitting in the energy region 120 MeV< Ui < 450 MeV has been 
performed recently, and the preliminary results are published in [18], [19]. Here we present 
the new data analysis for this experiment [20] . The results obtained confirm the existence 
of photon splitting phenomenon. They also make possible the quantitative comparison 
with the theoretical predictions. Moreover, the attained experimental accuracy allows one 
to distinguish between the theoretical predictions obtained with or without accounting for 
the Coulomb corrections. It turns out that the Coulomb corrections essentially improve 
the agreement between the theory and the experiment. 
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2. General discussion 



Let a photon with 4-momentum ki {k^ =^71 = 1^1!) and polarization vector ei turns in 
the electric field of an atom into two photons with 4-momenta ^2,3 (/C23 = ^2,3 = 1^2,31) 
and polarization vectors 62,3. We assume that the momentum transferred to a nucleus is 
small as compared to the mass of the nucleus. This condition allows us to neglect the recoil 
effects and consider an atom as a source of a time-independent electric field. Then, the 
final photons share the energy of the initial quantum: ui = UJ2 + UJ3. Later we will see that 
for arbitrary photon energy the main contribution to the total cross section comes from 
recoil momenta smaller than several tenth of electron mass. Therefore, the approximation 
of an external field is valid everywhere except the kinematic region where the differential 
cross section of the process is negligibly small. 

The most convenient way to take into consideration an external electromagnetic field 
in quantum electrodynamics is to use the Furry representation. In this approach the 
amplitude of a process is described by a set of Feynman diagrams where the electron 
lines correspond to the Green functions of the Dirac equation in the field. As a result, we 
obtain the amplitude in the form of series in the parameter a, where the coefficients are 
exact in the external field strength. 

Here we consider the amplitude M of photon splitting only in the lowest order in a which 
is given by the Feynman diagrams shown in Fig. 1. 




Fig. 1. Feynman diagrams for the photon splitting amplitude in the Furry representation. Thick 
Unes denote the electron propagator in an external field 

The corresponding analytical expression for M reads 



M = ie^ J ^ J dri dr2 dr^ exp[i{ki ■ vi - k2 ■ r2 - ■ r^)] (2.1) 
xTr[eiG(ri, r2\e - u;2)e;G'(r2, r3\e)elG{r3, ri\e + uj^)] + {k!^ ^ k^, e.2 ^ 63) . 

Here e = e^7^ = —67, 7^ being the Dirac matrices, G{ri,r2\e) is the Green function of 
the Dirac equation in the external electric field: 

G{ri,r2\e) = {ri\[V -m + t0]-^\r2) , P = 7°(£ - f/(r)) - 7 • p , (2.2) 



5 



p = — iV, U{r) is the potential energy of an electron. Note that, according to the Furry 
theorem, the amplitude (2.1) is an odd function of the external field strength (of the 
parameter Za for a Coulomb field). 



The differential cross section of photon splitting has the form 



da 



1 



Ml'^cul x{l - x) dx dn2 dfl, 



(2.3) 
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where x — uj2/oj\ (so that uj^ — ijJi{1 — x)), and ^2,3 are solid angles of ^2,3- If is convenient 
to perform the calculations in terms of the helicity amplitudes Mx^\^\^{ki,k2,kz)i with 
Aj = ±1. In this case the polarization vectors satisfy the relations e;^ ■ ^ = and e\ x 
k = i\uje\. In fact, it is sufficient to calculate three amplitudes, e.g., M+++(fci, ^2, ^3), 
M++_(fci, ^2, ^3), and M_| (fci, ^2, ^3)- Due to parity conservation and identity of pho- 
tons, other amplitudes can be obtained by the substitutions (see, e.g. [6]). The cross 
section da\^x^\^ for circularly polarized photons is given by (2.3) with M being Mx^x^x^- 
For unpolarized initial photon the cross section summed up over the polarizations of final 
photons is given by (2.3) substituting for |Mp 



Below we use the system of coordinates with 2;-axis directed along ki so that — a ■ 
ki/uji and a^ — a — ttzki/uji for an arbitrary vector a. Then the amplitudes depend on 
X = UJ2/UJ1, polar angles ^2,3 and azimuth angles 02, 3 of the vectors ^2,3. For spherically 
symmetric potential the quantity IMA1A2A3P depends on azimuth angles only via = 
4>2 — 4>3- In this case, due to the parity conservation, we have 



where Aj denote the hehcity opposite to Aj. Generally speaking, |MaiA2A3(0)P is not an 
even function of 0. 

Four theoretical approaches, having different ranges of applicability have been used for 
calculation of the cross section of photon splitting in an atomic field. Three of them use 
the first-order perturbation theory in the field strength, thereby implying Za <^ 1. They 
are 

• The approach exploiting the Heisenberg-Euler effective Lagrangian, which can be used 
only for cui <tC m. 

• The Weizsacker- Williams approximation, valid for ooi ^ m, and providing the logarith- 
mic accuracy for the cross section of the process. 

• The calculation of the amplitude exactly with respect to all kinematic parameters of 
the problem. This approach will be referred below as the Born approximation. 



(2.4) 



(2.5) 
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The foTirth method is based on the quasiclassical approximation vahd for oji ^ m. In 
contrast to the approaches itemized, it gives the cross section exact in the parameter Za, 
and provides a power accuracy in the small parameter m/cui. 

All these approaches will be reviewed below. 
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3. Low-energy photon splitting 



The Heisenber-Euler effective Lagrangian (HEL), derived in [21], describes self-action of 
an electromagnetic field due to the vacuum polarization when the momenta of quanta are 
small compared to the electron mass m. This Lagrangian depends on two invariants of 
the electromagnetic field: 

J^=-{H'^-E^) , g = E H. 
The lowest order term of the expansion of the HEL with respect to these invariants reads 



3607r2m4 



E' - Hy + i{E ■ ny 



(3.1) 



To calculate the photon splitting amplitude, it is necessary to present the electromagnetic 
field in (3.1) as a sum of a quantized field and the classical electric field of an atom and 
take the corresponding matrix element. The result obtained accounts for a single Coulomb 
exchange and holds for uji <^m (and, hence, uj2,z ^ fn) because the momentum transfer 
A = ^2 + ^3 — ki in this exchange is also small as compared to m. We emphasize that it 
is impossible to calculate the higher order corrections in Za (Coulomb corrections) to the 
photon splitting amplitude using the HEL. This is due to the fact that at the multiple 
Coulomb exchange the typical momenta of individual Coulomb quanta are of the order 
of m though the sum of these momenta A is small. 

The calculation of the low-energy photon splitting amplitude with the use of the HEL was 
performed in [7,8]. Unfortunately, these papers contain some misprints in the calculations 
of total cross section as was pointed out in [11], where the photon splitting was investigated 
numerically using the amplitudes obtained in [6] in the Born approximation. 

It follows from (3.1) that the amplitude of the process in an unscreened Coulomb field 
has the form 



M 



A1A2A3 



457r2m4A2 





ei)(e; 




■e;)(ei 


+(A 





l + A2A3 + -Ai(A2 + A3) 
7 

l-AiA3 + -A2(Ai-A3) 
7 

1 - A1A2 + -A3(Ai - A2) 



(3.2) 



In the case of the electric field of an atom this amplitude should be multiplied by an 
atomic form factor [1 — F(A^)] which accounts for screening. This form factor, being the 
Fourier transform of a charge density (in units of Z\e\), vanishes at A = and tends 
to unity for A —> 00. As a function of A it has a typical scale of = maz^/^. These 
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features are illustrated by a simple representation for the form factor suggested in [22] 

1-F(A^)=A^|:^^, (3.3) 

where 

= 0.1, 02 = 0.55, ag = 0.35, A = /?o&i, (3.4) 
6i = 6, 62 = 1.2, 63 = 0.3, /3o = 

For uoTc ^ 1, the effect of screening leads to the strong suppression of the cross section of 
photon splitting as compared to the case of unscreened Coulomb field. 

The amplitudes (3.2) obey the relation \Mx^\^x^\'^ — IMa^AjAsP (A denotes the hehcity 
opposite to A). Then, as follows from (2.5), the quantity is an even function of 

0. 

The angular distribution of the final photons is rather comphcated. This is illustrated in 
Fig. 2, where the cross section da / dx dD,2d^3 is shown as a function of 9^ at 62 — tt/IO 
and = 0, TT. 




0.5 1 1.5 2 2.5 3 



O3 

Fig. 2. The cross section da / dx d^l2dVL'i in units of K)-"^ Z'^a^ui^ /m^ = 3.09 • Z'^{iOi/mf pb as 
a function of 6z for x = 0.4, 62 = vr/10, and = 0, vr (dashed and solid curves, respectively). 

For (f) — TT there is a double peak with a narrow notch. Such a structure occurs in the 
range of the momentum transfer A <^ui where ^2,3 1. Then the bottom of the notch is 

at = TT, 63 = 62X/ (1 — x), when Aj^ = 0. Note that for = the differential cross section 
has a wide deep around = 92- In fact, such a suppression of the cross section, when the 
vectors ^2 and are almost parallel, persists for any ^2. The cross section da/dxdfl2 
differential with respect to the momentum of one photon is plotted in Fig. 3 for different 



9 



values of x as a function of cos ^2- It is seen that the photons with 1 — ,x ^ 1 (,t ^ 1) 
are emitted preferably along ki {—ki). It is also interesting to consider the dependence 
of the cross section on the azimuth angle (p = 4>2 — 4>3, shown in Fig. 4. This cross section 
is symmetric with respect to the replacement a; — > 1 — x. One can see that for all x the 
0-distribution has a minimum ai (p — 0. The cross section being the function of |0| has 
a maximum that moves from |0| = tt to |0| = 2.35 when x increases from to 0.5. The 
cross sections da \^\^\^ / dx for different helicities as well as da/dx for unpolarized photons 
are shown in Fig. 5. Note that da++_/dx is not symmetric with respect to substitution 
X — > 1 — x. Instead, after this substitution da^+^/dx turns into da^^^/dx. 




-0.75 -0.5 -0.25 0.25 0.5 0.75 1 
COS $2 

Fig. 3. The cross section da/dxdQ2 in the same units as in Fig. 2 versus cos02 for 
X = 0.2,0.4,0.6,0.8 (dotted, dashed, dash-dotted, and soUd curves, respectively). 
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Fig. 4. The cross section da/dx in the same units as in Fig. 2 versus for x = 0.1,0.3,0.5 
(dashed, dotted, and solid curves, respectively). 
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0.2 0.4 0.6 0.8 1 

X 

Fig. 5. The cross section da/dx in the same units as in Fig. 2, for unpolarized photons (sohd 

curve) and for special cases of heUcities: A1A2A3 = H (dotted curve), + H — (dashed curve), 

and H (dash-dotted curve). 

Numerical integration gives the following result for the total cross section in the Coulomb 
field 



rfa; = 7.6-10-^^ - , (3.5) 

2 J ax m"^ \m J 



which coincides with that obtained in [11]. The factor 1/2 in front of the integral accounts 
for the identity of photons. 
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4. Weizsacker- Williams approximation 



The main contribution to the total cross section of photon splitting for ui ^ m comes 
from the region of small angles ^2,3 ~ m/uji and can be obtained within the logarithmic 
accuracy using the Weizsacker- Williams (WW) approximation. Using this approximation, 
rather accurate estimate of the total cross section was first made in [9] and later in [7,8]. 
The partially integrated cross section as well as the total one were calculated in this 
approximation in [10]. In the WW approximation the interaction with a Coulomb quantum 
is replaced by that with a real photon having the 4-momentum q — —q^ki/ui). Then 
the cross section of photon splitting in WW approximation is expressed via the cross 
section of photon-photon scattering dcr^^ as follows 



Aeff 

da=^-^^L{s)da,,, L{s)= [^[1-F{A')]\ (4.1) 

S = {ki + qf = AUiq^ , Amin = — , Ae// = \/s, 

where [1 — F(A^)] is the atomic form factor (3.3). The factor in front of dajj{s) is the 
spectral distribution of equivalent photons. The photon-photon scattering cross section 
has the form 

dk dk 

da,, = \M\H{k^ + q-k2- ks) ^^ / , (4.2) 

where M is the amplitude of photon-photon scattering [2]. To eliminate the four-dimensional 
5-function in the cross section, it is necessary to integrate over the momentum of one of 
the final photons, and over qo in (4.1). We have k2± + = or = tt and ^262 = ^^3^3 
for small polar angles. Then the differential cross section da/dxdk2± has the form 

da _Z'a\M\'L{s) „ _ k^^ ^^^^^ 



dxdk2± 87r^a;(l — x)s^ ' x{l — x) 

The differential cross section for unpolarized photons is given by (4.3) substituting for 

\M 



2 



\Mf ^ ^[|M++++P + |M++-r + + |M+-+|2 + 4|M+++_|2] , (4.4) 

where Mx-^^x^x2\3 helicity amplitudes of photon- photon scattering which explicit 

form is given in Appendix A. 

Integrating (4.3) over all variables we obtain the total cross section of photon splitting 
which can be represented in the form [10] 
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2Z2a 7ds 



— J —(Tjj{s)L{s) . 



(4.5) 



Since the cross section ay^{s) falls off for s ^ w?, the integral in (4.5) converges at 
s w?. Within the accuracy provided by the WW method we can take the quantity 
L{s = m^) outside the integral. The magnitude of L{im?) strongly depends on the ratio 
of A~J„ = ujijw? to the radius of screening Tc ~ \moLZ^I'^\~^ . In the region ujx > m^rc 
we have L^m^) = In mrc = ln(Z~^/^/Q;) and a turns out to be independent of cui. For 
m <^ LUi < m'^rc, when the screening may be neglected, we have L{rn?) = ln(a;i/m). As a 
result, for the total cross section we have 

\n(uji/m) for m <^ a;i < mZ'^^^/a 

(4.6) 

ln(Z-V3/Q;) for ui > mZ-yya 



a ^0.725— ^ X < 



Recall that this result has the logarithmic accuracy. The coefficient 0.725 in (4.6) is taken 
from Eq. (4.2) of [11]. Instead of this number, the approximate value 0.753 was obtained 
in [10]. 

The detailed comparison of the results obtained within WW approximation, the Born 
approximation, and the quasiclassical approximation will be performed below. Here we 
only note that the accuracy of (4.6) at uji > mZ-i/V a is about 12% but for the differential 
cross section (4.3) this accuracy may be essentially worse. 
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5. Born approximation 



In two previous sections we liave considered photon splitting within the first-order pertur- 
bation theory in Za using some additional approximations valid in restricted kinematic 
regions. The analytical expression for the amplitudes in the Born approximation valid for 
arbitrary momenta of photons was obtained for the first time in [6] and later confirmed 
in [2]. The Feynman diagrams corresponding to the amplitude in the Born approximation 
are shown in Fig. 6. 





+ 



62^63 



Fig. 6. Feynman diagrams for the Born photon splitting amplitude. Wavy lines with crosses 
denote Coulomb quanta. 



For unpolarized photons the result of [6], multiplied by the atomic form factor, reads 

(5.1) 



dx dQ2 dQ, 



where 



C = -2k2h/m^ , p = 2hk2/m^ , 7 = 2kih/m^ , 7] = 4(fci ■ [fcs x fcg])'/^^ • (5-2) 

The quantities Xi and i/i depend on cuj and invariants (, (3, 7. Their explicit form being 
very cumbersome is presented in Appendix B. 

The numerical investigation of various differential cross sections of photon splitting, based 
on the analytical result of [6], was performed in [11]. The results of [11] are quoted in this 
section. Shown in Fig. 7 is the dependence of the total Born cross section on ui for 
screened and unscreened Coulomb potentials. The values of this cross section are listed 
in Table C.l of Appendix C for various Z and ui. 
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500 1000 1500 2000 2500 3000 

ui/m 

Fig. 7. The cross section (7/Z^(nb) plotted as a function of uoi for a Coulomb potential (solid 
curve), and the screened potential at Z = 1 and Z = 92 (dashed and dotted curves, respectively). 

The numerical results for the total cross section at uji < 0.1m, obtained in [11], coincide 
with (3.5). The analytical fit to the cross section in a Coulomb field, reproducing the cross 
section in the range 10m < coi < 3000m within 0.1%, reads 



22.372 In — - 47.313 + 82.536 31.64 — 

m a;i V a; 



nb . 



(5.3) 



In contrast to the case of a Coulomb field, when the total cross section grows logarith- 
mically with increasing cui, the photon splitting cross section in the case of a screened 
Coulomb field is saturated for ui ^ mZ"^!^ j a. This trend is clearly seen in Fig. 7. Note 
that the difference of the cross sections for screened and Coulomb fields at lo\ — mZ~^/^ /a 
is small (a few percent). 

The quasiclassical approach presented in the next Section allows one to obtain the photon 
splitting cross section for a; m with a power accuracy in the parameter m/oji. Therefore, 
the results obtained within this approximation reproduce for Za <^ 1 two first terms of 
(5.3). Thus, the accuracy of the quasiclassical approach can be estimated by the relative 
magnitude of two last terms in (5.3). This accuracy turns out to be better than 1% at 
> 70MeV. 

In Fig. 8 the spectrum a~^da/dx for a Coulomb field is plotted as a function of x for 
different values of (Vi. One can see that the shape of the spectrum changes significantly 
with increasing photon energy. The detailed data obtained in [11] for spectrum in screened 
and Coulomb fields are presented in Appendix C in the Table C.2. 

In order to realize the accuracy of the Born approximation at Za ~ 1 it is necessary to 
calculate the Coulomb corrections. This issue is elucidated in the next Section. 
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0.2 0.4 0.6 0.8 1 

X 

Fig. 8. The spectrum a^^da/dx is plotted for the case of a Coulomb potential as a function of 
X for wi/m = 1, 10, 100 (solid, dashed, and dotted curves, respectively). 
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6. High-energy photon spUtting 



Small scattering angles of particles provide the applicability of the quasiclassical approach 
to the description of various high-energy processes in an external field. This approach 
gives the transparent physical picture of phenomena and essentially simplifies calculations. 
Recall that the region of small angles ^2,3 ^ 1 makes the main contribution to the total 
cross section of photon splitting for 3> m. That is why we confine ourselves on this 
kinematic region when considering the process in this Section. Below we quote the results 
obtained in the quasiclassical approximation in [13-15] . The technique used in these papers 
was developed in [23,24] for the investigation of coherent photon scattering in a Coulomb 



We choose the origin at the center of a potential and direct the 2;-axis along ki. Then 
characteristic values of transverse components of Vi making the main contribution to the 
integral in (2.1) at given momentum transfer A = ^2 + ^3 — ki arc of the order of 
p ~ 1/A. At small angles the characteristic angular momentum is / ~ a;/A ^ 1, and the 
quasiclassical approximation can be applied. 

It is important that for uj ^ m the exact in Za amplitude of photon splitting in the 
atomic field can be obtained from that calculated for a Coulomb field. For this purpose, 
wc represent the amplitude as a sum of the lowest in Za term (Born amplitude) and 
the Coulomb corrections. As explained above, the Born amplitude of the process in the 
atomic field is a product of that in the Coulomb field and the atomic form factor. In the 
perturbation theory the Coulomb corrections correspond to the multiple interaction with 
quanta of an external field. The sum of the momenta of these quanta equals A. It turns 
out that even for A <^ m the characteristic magnitude of the momentum of each quantum 
is not small as compared to m. As m 3> all form factors appearing in the calculation 
of the Coulomb corrections should be set to unity. Thereby, the Coulomb corrections are 
the same for the atomic and Coulomb fields. 



6. 1 Quasiclassical picture of the process 

The calculation of the amplitude of photon splitting is essentially simplified by using the 
Green function D{ri,r2\s) of the "squared" Dirac equation: 



The Green function G{ri,r2\e) of the Dirac equation can be obtained from D{ri,r2\s) 
by means of the relation 



field. 



D{ri,r2\e) = {ri\[V^ 



(6.1) 



G{ri,r2\e) = [7°(£ - [/(n)) + ^7 • Vi + m]D(ri, raje) . 



(6.2) 
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As shown in [13,14], the initial expression for the photon sphtting amplitude (2.1) can be 
presented as a sum of two contributions, containing either three or two Green functions 
D{ri,r2\s): M = M^^) + M^^). The term M^^) is given by 



^(3) = -e^ y — j dridr2dr3 exp[i(fci • ri - ^2 • r2 - fcs • rs)] 

xTr{[(-ei^i - 2ei ■ p)D(ri, r^le - cu2mefk2 - 2el ■ p)D{r2, r^\e)] 
)<[{elk2.-2el- p)D{r^,ri\e + uj^)\] + (fc^ ^ /c^,e2 ^ 63). 



(6.3) 



Here the operator p = — iV differentiates the Green functions D with respect to their 
first argument. The term M^'^^ reads 



M(2) ^ie^ J dridr2 TV{exp[i(fci • n - ^2 • r2 - feg • r2)] 

xe^ • e3[(-eiA;i - 2ei ■p)D{ri,r2\e - t^i)]D(r2, nle:) 



+ 



exp[i(fci ■ ri - ^2 ■ rs - fca ■ ri)]ei ■ e3D(ri, r2|e - UJ2) 



x[(e;^2 - 2e; • p)D{v2, Ty\e)\ + (A;^ ^ A;^, 62 ^ 63)] }. 



(6.4) 



Below the Green function D(rx^r2\s) will be called the "electron Green function" for 
£ > and the "positron Green function" for e < 0. According to [23,24], for high energies 
the main contribution to the integrals in (6.3) and (6.4) comes from the region of inte- 
gration over the variables (^^-coordinates of r^) in which z' < z for the electron Green 
function D(r, r'\e) and z' > z for the positron Green function. In terms of non-covariant 
perturbation theory, this corresponds to the contribution of intermediate states (see Fig. 
1) for which the difference between their energy and the energy of the initial state 
Eq = uji is small as compared to £'0. 

The lifetime of the intermediate state in photon splitting can be estimated with the help 
of the uncertainty relation as r ~ uji/{m^ + k\)^ where k_[_ = max(|fc2_L|, 1^3x1)- The char- 
acteristic transverse distance between the virtual particles (transverse size of the loop) 
can be estimated as (m^ -|- kj_)~^^'^, which for uji ^ max(m, k±) is much smaller than the 
length of the loop r. The quantity p ~ 1/A introduced above can be interpreted as a typi- 
cal transverse distance between the loop and the origin. In the small-angle approximation 
(^2,3 -C 1) we have for 

= Ai + = (fc2x + fesx r + \ f ^ + • (6.5) 

4 V u;2 a;3 y 



It is seen from (6.5) that Aj^ 3> A^ everywhere except a narrow region where Aj^ -C k^. 
As pointed out in [14], the amplitude at A^ ~ A;^ can be obtained from that derived for 
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A± A2 by a simple substitution that will be discussed below. The condition ^ 
implies that p/r ^ 1 and, therefore, the angles between fcj and in (6.3) and (6.4) 
are either small or close to vr. The additional restrictions on the region, making the 
main contribution to the amplitude for A^ ^ A^,, follow from the explicit form of the 
quasiclassical Green function. For the term M^^-* (6.3) they are 

zi < Z2,Z3 , Zi < , max(2;2, ^3) > 0. 



Similarly, the main contribution to the term M*-^-* (6.4) is given by the region zi < and 
Z2 > 0- All these conditions allow one to depict the main contribution to the amplitude 
M^^^ and M*^^^ in the form of diagrams, shown in Fig. 9. 




Fig. 9. Diagrams corresponding to the terms M^^^ (left column) and M^'^^ (right column). 

The explicit form of vertices is obvious from (6.3) and (6.4). The electron Green functions 
are marked with left-to-right arrows, and positron ones with right-to-left arrows. The 
arrangement of the diagram vertices is space ordered. With the use of these diagrams one 
can easily determine the limits of integration over the energy and coordinates. Diagrams in 
Fig. 9 have a transparent interpretation. For example, the upper left diagram corresponds 
to the following picture: the photon with momentum fci produces, at the point ri, a pair 
of virtual particles which is transformed at the point r2 into a photon with momentum 
^2. Between these two events the electron emits a photon with the momentum ^3 at the 
point r^. 

The quasiclassical Green function G{ri, r2\e) of the Dirac equation in a Coulomb potential 
was first derived in [23] starting from the exact Green function of the Dirac equation 
[25]. For the case of almost collinear vectors Vi and r2 the quasiclassical Green function 
in arbitrary spherically symmetric decreasing potential was found in [24]. Later in [26] 
the quasiclassical Green function was obtained in an arbitrary localized potential which 
generally possesses no spherical symmetry. 

As explained above for the calculation of the amplitude of photon splitting in atomic 
field it is sufficient to know the quasiclassical Green function D{ri, r2\e) of the "squared" 
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Dirac equation in a Coulomb field for the case of almost coUinear vectors fi and r2- The 
corresponding expression derived in [13,14] has the form 



D{ri,r2\s) 



le 



m(ri+r2) 



87r^Krir2 



/ dq 


a ■ a 







X exp 



2 £ rir2 



q2 



iZarj 



(6.6) 



for a small angle between vectors r-i and —r2- Here ot — 7^7, — — w? , rf — sgn(£), 
''^1,2 — (»'i,2)±/'"i,2, and qr is a two-dimensional vector lying in xy plane. The direction 
of z axis is arbitrary provided that ?92.3 ^ 1- Expression (6.6) contains only elementary 
functions, and the angles t?i and t?2 appear only in the factor exp[iqr ■ (t?i + 'i?2)]- For a 
small angle between vectors ri and r2 we have for D{ri, \ e) 



D{ri,r2\ e) 



giK|ri— r2 



47r|ri - r2| \r2J 



\ iZar]sgn{ri—r2) 



(6.7) 



One can see that the Green function (6.7) differs from that at vanishing potential only by 
a phase factor. It is easy to check that after the substitution of the expressions (6.6) and 
(6.7) for the Green functions to the splitting amplitudes (6.3) and (6.4) all phase factors 
of the form r^*^° cancel. 



6.2 Amplitude of the process 



At the calculation of the amplitude it is convenient to eliminate the z component of the 
polarization vectors Cj using the relation Ci ■ ki — which leads to e^^ = — e^j, • ki±/kiz. 
After that, in the small- angle approximation we can neglect the difference between the 
polarization vector e2± (esj,) and that of photon, propagating along the z axis and having 
the same helicity. As a result, we obtain the amplitudes M\-^^x,^x,^{ki, k2, k^) containing 
only the transverse polarization vectors e and e* for positive and negative helicities, 
respectively. Though the Green function D{ri,r2\e) (6.6), (6.7) has a simple form, the 
integration in (6.3) and (6.4) was very sophisticated. The final expressions, obtained in 
[15], for the helicity amplitudes exact in the parameter Za have the form 
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u;2^^3 
e G 



v. 



[4 + EK^) B + 2k3 [4 + (e* • 6>23) (e • b) 

a;2«3-B 



6(;2G'i e ■ 6 a;2Gi e • 6 e - G 



+m {u;,u;2 - 2«2«3)j , , , ^^^^ 

+2^3 (k^ + (e* • 62) (e • 6) + {UJ2UJ3 - 2£K3) 

where the following notation is used 

K2 = ti^2 - £ , K3 = ^3 + £ , O23 ^ O2 — 03 , + , B ^ + b'^ , 

a^X- + K202, b^X+-K303, C ^ X+ - 8023 , X±^X±A±/2, 

N — 8e^ Z a[7r^ A^UJi0i)2UJ3]~^ , OjIVi — (k,2X+ + K3X^Y ~ <^2<^3'«2'«3^23 ~ ^0 ' 

a;|X>2 = {^^+ ~ K'2,X-)'^ — 0JiUJ2K'3e9\ , oulT>3 — {k2X+ + £a;_)^ + CU1CU3K2S9I . 
The functions G and Gi are 
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PlzM , (6.9) 



where c± = m^ + a;^, ^ = 2m^A^/(c+c_) — 1, and P/^alO = (iPiZa(0/'^C is the derivative 
of the Legendre function. The amplitudes with other hehcities can be obtained from (6.8) 
using the following relations 



M+_+(fci, fc2, ^3) = M++_(fci, fcg, ^2), 
M_A2A3(fcl,fc2,fc3) =M+A2A3(fel,fc2, ^3) (c ^ 6*), 



(6.10) 
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where A denotes the hehcity opposite to A. 



Since the functions G and Gi are independent of the energy e, the integrands in (6.8) 
are rational functions of £, in which all the denominators are quadratic forms of this vari- 
able. Therefore, the integrals over e can be expressed via elementary functions. Resulting 
formulae being rather cumbersome are not presented here explicitly. Performing the inte- 
gration over e in (6.8), one obtains a twofold integral over x for the amplitudes of photon 
splitting. It has been checked numerically that in the limit Za our results (6.8) agree 
with those, obtained in [11] in the Born approximation. 

Let us comment on the appearance of the total momentum transfer squared (see 
(6.5)) in the coefficient in (6.8). The amplitudes (6.8) were derived at the assumption 
Aj^ » A^ when there is no difference between A^ and A^^. Nevertheless, if we retain the 
term A^ in A^ in the coefficient A", the formulas (6.8) become valid also for Aj^ ~ A^. 
This can be explained as follows. In the region A^ ~ A^ the Born amplitude prevails 
over the Coulomb corrections (cf. (D.4) and (D.5)). On the other hand, in the small-angle 
approximation it contains the momentum transfer in the form A^/A^, where A is just 
the total momentum transfer. 

For 0^2^25 i^3^3 ^ we can neglect the electron mass in the expressions (6.8). In this case 
(zero mass limit) the amplitudes are simplified, which allows us to perform the further 
analytical investigation of the process (see Appendix D). 



6.3 Cross section 



In the small-angle approximation (^2,3 ^ 1) the cross section (2.3) can be represented in 
the form 

dxdk,^dk,^ (6.11) 
' ' 2%5^fx(l-x) ^ ' 



In this subsection a screened Coulomb potential is used in the numerical calculations. 
As was explained above, to take screening into account, the Born part of the amplitude 
should be multiphed by an atomic form factor, e.g. (3.3). To illustrate a magnitude of the 
Coulomb corrections, the exact and Born differential cross sections for unpolarized photons 
da/dxdk2j_dk3± arc plotted in Fig. 10 as a function of k2±/m for A;3_l = 2m, x = 0.1, 
cui = IGeV and the azimuth angle (the angle between the vectors k2± and ^3^) = 
and TT. The value Z = 83 (bismuth) was chosen since bismuth atoms determine the cross 
section of photon splitting in the experiment [18]. A wide peak for azimuth angle (f) — ir 
is due to small momentum transfer A. There is a narrow notch (cf. Fig. 2) in the middle 
of this peak ( at k2±/m = 2 ) where the condition A_l = k2± + k^^ = is fulfilled. Such a 
structure was discussed first in [12]. The width of the notch is about max(A^/m, /9o/m). 
Recall that A^ is the longitudinal component of the momentum transfer defined by (6.5), 
and Po — mZ^/^/121 (3.3) characterizes the effect of screening. In our example Po is 
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Fig. 10. Differential cross section da / dxdk2j_dk3± vs k2±/m in a screened Coulomb potential for 
different azimuth angle <f> between vectors k2± and fcs^; ^ = 83, x = 0.1, uii = IGeV, k^^ = 2m. 
The dashed curve (Born approximation) and the solid curve (exact cross section) correspond to 
(j) = IT. The dash-dotted curve (Born approximation) and the dotted curve (exact cross section) 
correspond to (j) = 0. 

larger than , so the width of the notch is roughly /3o/m = 3.6 x 10~^. Let us note 
that for uJi ^ m the differential cross section, expressed in terms of k2±, k^^, x, and Ui, 
depends on uJi only via A^. Due to this fact, the differential cross section is independent 
of Ui in the region where we can neglect A^, i.e. outside the notch vicinity. For small 
A as compared to k± the Born amplitude, which is proportional to 1/A, is much larger 
than the Coulomb corrections. That is why the exact and Born cross sections coincide 
within the peak region. Outside this region the Coulomb corrections essentially modify 
the cross section. The points of discontinuous slope of the curves in Fig. 10 are related to 
the threshold conditions for the production of real electron-positron pair by two photons 
with the momenta k2 and k^: 



It is worth noting that the cross section da\^x^\^ / dx dVL2dVL^ for definite helicities of pho- 
tons is not an even function of the azimuth angle 0. Instead, due to the parity conservation 
the following relation holds 



As an illustration of the azimuth asymmetry of the cross section, the dependence of 
and ^ = 2(5*+ — S'_)/(S'+ + SJ) on the azimuth angle is shown in Fig. 11, where S± = 
SA2A3 da±x^x^/dx dQ2d^3 are the differential cross sections for positive (5*+) and negative 






(6.13) 
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{S-) hclicities of the initial photon. Again, the exact in Za cross section (sohd curve 
in Fig. 11a) is much smaller than the cross section obtained in the Born approximation 
(dashed curve ibid). The magnitude of azimuth asymmetry (Fig. 116) reaches tens of 
percent. For the parameters used in our example the effect of screening is not essential. 




Fig. 11. The dependence of S-^- (a) and ^ = 2(S^ — 5'_)/(5+ + S-) (5) on the azimuth angle (p at 
Z = 83, LOi = IGeV, x = 0.3, 02 = 0.014, and 6*3 = 0.02. Here S± = J^x^x-, da±x2Xi/dxd^2d^2, 
are the differential cross sections for positive (5*+) and negative {S-) helicities of the initial 
photon. Solid curves represent the exact in Za results, dashed curves are obtained in the Born 
approximation. 

In Fig. 12 the differential cross section m^aQ^da/dx dk2± is shown depending on k2±/m 
for a screened Coulomb potential at coi/m — 1000, Z — 83, and different x, 

gg^ ^'^ff ^ 0.782 -10-^ Z^b. 

Solid curves represent the exact cross sections, the dashed curves give the Born results, 
and the dotted curves correspond to the WW approximation (4.3). The cross section 
exhibits a thresholdlike behavior in the vicinity of the point k2± — kfh — 2y^x(l — x)m, 
where both conditions Aj_ — k2± + ks± — and (6.12) hold. Under these conditions 
the peak in the cross section da / dxdk2i_dk^±_ seats on the boundary of the kinematic 
region in which the production of a real electron-positron pair by two photons with the 
momenta ^2 and ^3 is possible. The cross section da/dxdk2± drops rapidly for ^ m 
(cx l/^lx)- As s^6^ Fig- 12, the Coulomb corrections to the cross section integrated 
over k3± noticeably diminish the magnitude of the cross section. Above the threshold 
{k2± > kth) the Coulomb corrections reach tens of percent while below the threshold 
the exact cross section is several times smaller than the Born one. It can be explained as 
follows. Above the threshold the main contribution to the cross section da/dx dk2± is given 
by the integration region with respect to fe3_L where max(A;|_|_/a;i, /3o) <^ A <^ k2±- As a 
result, the Born cross section is logarithmically amplified as compared to the Coulomb 
corrections. Far below the threshold where k2± <C kth, the region A <g k2± does not 
contribute to the cross section d(T/dxdk2± since in this region A;3_|_ ^ k2_L ^ m, and the 
amplitude is suppressed as a power of A;|_|_/m^. Therefore, below the threshold the main 
contribution to the cross section is given by the region k3± ~ m, where the amplitude, 
exact in Za, drastically differs from the Born one. Due to the same reason below the 
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Fig. 12. m^fjQ ^ da / dxdk2± vs k2±/m for a screened Coulomb potential, LO\jm = 1000, Z = 83, 
£70 = a^{Zaf/{A'K'^w?) = 0.782 • lO'^ Z^b, x = 0.2 (a), 0.5 (6), and 0.8 (c). The dashed curves 
correspond to the Born approximation, the solid curves give the exact result, and the dotted 
curves show the results of the WW approximation. 

threshold the cross section in the WW approximation (dotted curves) differs essentially 
from the Born result because the former comes from the region A -C /c2±. As follows 
from (A. 3), for k± m the amplitudes entering the cross section (4.3) are suppressed 
as k\/m^. One can see from Fig. 12 that in accordance with the expression for ktu the 
position of the peak is the same for a; = 0.2 and x = 0.8. Nevertheless, the magnitude of 
these two cross sections are significantly different, especially below the threshold. However, 
as expected and verified numerically, the cross section da/dx at x = 0.2 coincide with 
that at X = 0.8. 

Let us consider now the magnitude of the Coulomb corrections to the cross section da/dx 
integrated over the transverse momenta of both final photons. The main contribution 
to this cross section is given by the region where A;2±i ~ n^- The Born contribution 
to da/dx contains large logarithm resulting from the integration over small momentum 
transfer region max(/5o, m^/a;i) <g A <^ m. For /3o ^ w?/uji, the cross section da/dx is 
independent of c^i, while for (5q <^ rr? /ujx it grows slowly (as ln(a;i/m)) when a;i increases. 
Since the Coulomb corrections to da/dx are determined by the region of momentum 
transfer A m, they do not depend on uji for uji ^ m . They also arc insensitive to the 
effect of screening. In Fig. 13 the exact (solid curve) and the Born (dashed curve) cross 
sections a^^da/dx are plotted as functions of x at cui/m — 1000, and Z — 83. As should 
be, the curves are symmetric with respect to the replacement a; — > 1 — x. The dotted curve 
shows the Coulomb corrections (the difference between the exact in Za cross section and 
the Born one) taken with the opposite sign. Note that their dependence on x is rather 
weak. The cross section da/dx in Fig. 13 increases rapidly when x approaches the points 
X — and x — 1. However, this cross section should vanish at x — and x — 1 due to 
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the gauge invariance of QED. It turns out that the decrease of the cross section da/dx 
occurs in the very narrow interval 5x ~ m} lio\ close to the points x = 0, 1. Therefore, the 
contribution from these domains to the total cross section a is negligible. In our example ( 
Z = 83, Ui/m = 1000 ), the exact result for a is 3.9 x 10~^b while the Born approximation 
gives 4.8 X 10~^b, the difference being 23%. 
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Fig. 13. The dependence of a^^ da/dx on x for a screened Coulomb potential, ivi/m = 1000, 
Z = 83. The dashed curve corresponds to the Born approximation, the solid one gives the exact 
result, and the dotted curve shows the difference between the Born cross section and the exact 
one. 

In Fig. 14 the Coulomb corrections to the cross section daddx divided by ctq = 0.782 • 
10~^Z^b are shown as a function of Z for x = 0.7 (solid line). Since their dependence on 
X is rather weak ( see Fig. 13), this curve allows one to estimate the magnitude of the 
Coulomb corrections for any x. To elucidate the role of higher-order Coulomb corrections 
we plot also the first Coulomb correction to the cross section divided by ctq (dashed 
curve). This ratio is proportional to Z^. As seen in Fig. 13, for sufficiently large Z the 
first Coulomb correction is irrelevant to the complete result for these corrections. 

To conclude this Section, we emphasize that the proper description of photon splitting 
in the electric field of a heavy atom requires the exact account of this field. For large Z, 
the contribution of the Coulomb corrections is always noticeable, though the magnitude 
of the effect depends on the type of a cross section and kinematic conditions. 
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Fig. 14. The dependence of the Coulomb corrections UQ^duc/dx on Z for a; = 0.7; sohd curve: 
complete result, dashed curve: first correction. 
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7. Experimental investigations of photon splitting 

The observation of photon spHtting is a difficult problem due to a smallness of its cross 
section as compared to those of other processes initiated by the initial photons in a target. 
The following background processes are significant: double Compton effect off the atomic 
electrons (je 77e), and the emission of two hard photons from e+e~ pair produced 
by the initial photon. The relative importance of these processes depends on the photon 
energy. For the energy uji ~ m, where the search of photon splitting was undertaken 
in two experiments [16,17], only double Compton scattering determines the background 
conditions. In these experiments, the photons from an intense radioactive source {Zn^^ 
with ui ~ 1.1 MeV in [16], and Co^° with cji ~ 1.17, 1.33 MeV in [17]) were used. The 
combination of the coincidence and energy-summing detection technique was applied. The 
number of events considered as candidates for photon splitting exceeded the theoretical 
expectations by the factor of 300 in [17], and by the factor of 6 in [16]. The role of the 
background process in [16] was unclear. Nevertheless, the authors of [16] concluded that 
their " results for Z ~ 28 strongly indicate the existence of photon splitting" . 

At high photon energy uj ^ m the emission of hard photons from e^e~ pair becomes 
most important as a background process. In 1973 the experiment devoted to the study 
of Delbriick scattering of photons in energy region 1 7 GeV was performed [28]. The 
bremsstrahlung non-tagged photon beam was used. Some events were assigned by authors 
of [28] to the photon splitting process. As shown in [10,29], the theoretical value for the 
number of photon splitting events under the conditions of the experiment was two orders 
of magnitude smaller than the experimental result. It was also argued that the events 
observed could be explained by the production of e^e~ pair and one hard photon. 

The first successful observation of photon splitting was performed in 1995-96 using the 
tagged photon beam of the energy 120 -^ 450 MeV at the VEPP-4M e^e~ collider [30] in 
the Budker Institute of Nuclear Physics (Novosibirsk). Another goal of this experiment 
was a study of Delbriick scattering [5]. The total statistics collected was 1.6 ■ 10^ incoming 
photons with BGO (bismuth germanate) target and 4 ■ 10* photons without target for 
background measurements. The preliminary results were published in [18], [19]. Here we 
present the last data analysis for this experiment [20]. 

The experimental setup is shown in Fig. 15. Some ideas of this setup were suggested 
in [31]. The main features of the experimental approach are: 

• The use of high-quality tagged photon beam produced by backward Compton scattering 
of laser light off high-energy electron beam. Thereby, the energy of the initial photon 
is accurately determined. 

• Strong suppression of the background processes by means of the detection of charged 
particles produced in the target and in other elements of the photon-beam line. 

• The detection of both final photons to discriminate the photon splitting events from 
those with one final photon produced in Compton or Delbriick scattering. 
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• The requirement of the balance between the sum of the energies of the final photons 
and the energy of the tagged initial photon. This provides the additional suppression 
of the events with charged particles missed by the detection system. 




Fig. 15. Principal scheme of setup: LKr calorimeter (1); scintillation veto-counter (2); beam 
dump (BGO) (3); He-filled tube (4m length) (4); target (BGO) (5); active collimator (BGO) 
(6); lead absorber (7); guiding tube for the gamma-quanta beam (8); cleaning magnet (9); passive 
lead collimator (10); TS scintillattion counter (11); Nd:YAG is the laser; TS1-TS4 are tagging 
system hodoscopes; Ml and M2 are bending magnets; LI and L2 are quadrupole lenses. 

At high energy of the initial photon uoi ^ m, the photon splitting cross section is peaked 
at small angles between momenta of all photons (~ m/uji). Therefore, a good coUimation 
of the primary photon beam is required. The ROKK-IM facihty [32] is used as the in- 
tense source of the tagged 7-quanta. The electron energy loss in the process of Compton 
scattering of laser light is measured by the tagging system (TS) [33] of the KEDR detec- 
tor [34]. The TS consists of the focussing spectrometer formed by accelerator quadrupole 
lenses and bending magnets, and 4 hodoscopes of the drift tubes. High-energy photons 
move in a narrow cone around the electron beam direction. The angular spread is of the 
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order I/7, where 7 = E^jeam/fn is the relativistic factor of the electron beam. The photon 
energy spectrum has a sharp edge at 



47 ^laser 



(7.1) 



1 + A-fUJiaser/m ' 



that allows one to perform the absolute calibration of the tagging system in a wide energy 
range. In experiment the laser photon energy was ujiaser — 1.165 eV, the electron beam 
energy Eheam = 5.25 GeV, and Ufh = 450 MeV. The photon energy resolution provided by 
the tagging system depends on the photon energy and on the position of scattered electron 
in TS hodoscope: it was 0.8 % at a;i = 450 MeV (at the center of the hodoscope) and ~ 5 % 
at cji — 120 MeV (at the edge of the hodoscope). The coUimation of the photon beam 
is provided by two collimators spaced at 13.5 m. The last collimator, intended to strip 
off the beam halo produced on the first one, is made of four BGO (bismuth germanate) 
crystals as shown in the separate view in Fig. 15. After passing through the coUimation 
system, the photon beam hits IXq thick BGO crystal target. 

Since one has to separate incident photons from those scattered in the target, certain 
angular region around the photon beam direction < is enclosed by the dump made of 
13 Xq thick BGO crystal installed in front of the photon detector. The incoming photons 
which passed the target without interaction are absorbed by the dump, and the only 
photons to be detected are those scattered outside the dump shadow. For too large value 
of many photon splitting events would be lost while for too small value of the 
non-interacted photons would overwhelm the data acquisition system. As an optimum, 
the value 9o = 2.4 mrad was chosen. 

All active elements used in beam line (collimators, target, dump, scintillating veto counter) 
set a veto signal in the trigger and their signals are used in the analysis for background 
suppression. The information from the target and beam dump is also used for measurement 
of the incoming photon flux. The liquid-krypton ionization calorimeter is used for the 
detection of the final photons. Its three-layer double-sided electrode structure enables 
one to get both (X and Y) coordinates for detected photons. The energy resolution of 



the calorimeter is 2%/ ^^Juj{Ge\f). The liquid-krypton calorimeter is described in details 
in [35,36]. 

In the experiment, the detected final photons had the polar angles in the region 2.4 mrad < 
^ < 20 mrad. The corresponding cross section is called "visible". Fig. 16 shows the 

calculated energy dependence of the total (a) and visible (6) cross sections for various 
processes initiated by photons in BGO target: e+e^ pair production, Compton scattering 
on atomic electrons, Delbriick scattering, and photon splitting. The calculation of the 
photon splitting cross section was performed using the results obtained in [13-15]. 

In the event selection procedure the following constraints were applied: 

• The absence of the signal caused by charge particles in all active elements of the photon- 
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Fig. 16. The calculated energy dependence of the total (a) and the visible (b) cross sections 
of various processes initiated by photons in BGO target (in units of barn per one molecule 
of BuGe^Oi2). 

beam line. 

• The balance of the tagged initial photon energy and the energy measured in the 
calorimeter within 3cr of its energy resolution. 

• The existence of two separate tracks at least for one (X or Y) coordinate in the calorime- 
ter strip structure. The track is found if there are close clusters in different layers. 

The fulfillment of the latter requirement strongly suppresses the contribution of the pro- 
cesses with one photon in the final state which could imitate two photon events in the 
calorimeter. 



The typical event which meets selection criteria is shown in Fig. 17. In this example two 
tracks are seen in both X and Y directions. The conversion of the first photon occurs in 
the Layer 1 while the second photon converts in the Layer 2. 
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Fig. 17. Energy profile in the calorimeter strip structure (left) and reconstructed kinematics 
(right) for a typical candidate to the photon splitting event. 

The experimental results are presented in Tab. 1 and in Figs. 18, 19 together with the re- 
sults of Monte-Carlo simulation based on the exact photon sphtting cross section. The en- 
ergy spectrum of the initial photons measured in the tagging system is shown in Fig. 18(a). 
Tab. 1 and Fig. 19 present the data summed up over the initial photon spectrum. The 
errors shown in the Tab. 1 arc statistical ones. The systematic error is determined by the 
accuracy of the measurement of the number of initial photons and by the uncertainty in 
the reconstruction efficiency of photon splitting events. The estimation of these systematic 
errors gives 2 % and 5 %, respectively. 

Table 1 

The number of reconstructed events. Here Q is the number of incoming photons. The quantity 
-^</><i50° is the number of events with the complanarity angle (j) < 150° (sec Fig. 19), iV^>i5o° is 
the number of events with cp > 150°. The quantities A^0<i5o° and N^^i^qo are normalized to the 
experimental statistics collected with the target. MC means Monte-Carlo simulation. 



DATA 


TARGET 


Q, lO'^ 


N<)i,>i50° 


N0<15O° 


Experiment 


Bi4Ge30i2 


1.63 


336±18 


82±9 


Experiment 


no target 


0.37 


10±3 


10±3 


MC photon splitting 


Bi4Ge30i2 


6.52 


364±10 


72±5 


MC Delbriick scattering 


Bi4Gc30i2 


1.63 


2±1 


16±4 


MC other backgrounds 


Bi4Ge30i2 


1.63 





16±4 



As seen from the Table 1 and Fig. 19a, the main part of the photon splitting events 

has, in agreement with the theory, a complanarity angle (p (the azimuth angle between 
final photon momenta) close to 180°. The choice of the interval (f) > 150° allows us to 
improve the signal-to-background ratio (see, e.g., the last two rows of the Table 1). Just 
this 0-interval was used to plot the distributions over polar angles in Fig. 19 and the 
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dependence of the number of reconstructed photon sphtting events on the initial photon 
energy E^s in Fig- 18(b). Note that for most of the events in this 0-interval, the variable 
^ = &min/ {djnin + ^maa;) is approximately equal to the ratio min(u;2, ^z) l^i since the main 
contribution to the cross section comes from the region \k21_ + k^i] -C /c2±, k^^i i-e. 
180° and a;2^2 ~ oj^9^. 



I 200 - 

150 - 

100 r- 



150 200 250 300 350 400 450 

Q Ets , MeV 




Fig. 18. (a) The photon energy spectrum measured in the tagging system (TS). (b) The number 
of reconstructed photon splitting events as a function of the tagged photon energy Ets- In 
plot (b) black circles present the experimental results, histogram is the result of Monte-Carlo 
simulation based on the exact in Za photon splitting cross section. 

The results presented in the Table 1 and in Figs. 18(b), 19 show a good agreement between 
the theory and the experiment. More precisely, the total number of reconstructed events 
in the experiment (see Table 1) differs from the result of MC simulation by 1.5 standard 
deviations. 



In order to demonstrate the role of the Coulomb corrections under the experimental 
conditions, we show in Fig. 20 the visible cross sections calculated exactly in Za and 
in the Born approximation as a function of the initial photon energy. For all energies 
considered, the Born result exceeds the exact one by 20 % approximately. Therefore, the 
use of the Born cross section for MC simulation would lead to the disagreement between 
the theory and the experiment of 3.5 standard deviations. In other words, the experimental 
results are significantly closer to predictions of the exact theory than to those obtained 
in the Born approximation. 

We conclude that the high-energy photon splitting in the atomic fields is reliably observed 
and well investigated. This nonlinear QED process is also studied in detail theoretically. 
At present, the experiment and the theory are consistent within the achieved experimental 
accuracy. 
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C 0^,, , mrod d 0mIn/(Qmin+0mcx) 

Fig. 19. The number of the selected photon sphtting events as a function of 
the azimuth angle between momenta of the outgoing photons (a); 
the polar angle 9min = mm{92,93} (b); 
the polar angle 9max = max{6'2,6'3} (c); 
the variable x = 9min/{9min + 9 max) (d). 

In figures (b), (c), and (d) only events satisfying the complanarity cut ^ > 150° (see plot 
(a)) are included. Black circles present the experimental results, histograms are the results of 
Monte-Carlo simulation based on the exact in Za photon splitting cross section. 
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Fig. 20. The visible photon splitting cross section calculated exactly in Za (1) and in the Born 
approximation(2) as a function of the initial photon energy. 
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Appendix A 



In this Appendix we present the result of [1] for the amphtudes of photon-photon scat- 
tering entering (4.4) 



C 

2(^^+f) 8 

e C 



2 + ^)g(7) 



™+r(7)]-^(l + ^j/(C,/3) 



1 



2(3- + ^^) 16 4 4 8' 

^ — H \ \ 

C 7 tu 



M+++_ = 8a^|l - 4 + - + -J [T(C) + T(/3) + T{^)] 
M,,__ = 8«^{l-i7(C,/5)-A,(C,,)-l/(/^,,)j, 



(A.l) 



Here 



g(.) = Jl7-3_arcsinh(^^ 



— 1 , = arcsinh^ f ^*^ ^ - 



c 





2/^2^3 



1 /• . ln[l + a;r(l - r) - iO] + ln[l + yr(l - r) - lO] 



m?x(l — a;) 



t(1 -r) {x + y)/xy 
_ 2kik2 _ kl^ 

1 P 



7 



2kik^ 



A-2 



m^fl — x) 



(A.2) 



In the expressions for C, and 7 we took into account the smallness of ^^2,3 and the relation 
feax = — fc2±, which is valid in the region of applicability of the WW approximation. 

The asymptotics of (A.l) at |CL/5)7 ^ 1 reads 



++++ ■ 



lla^ 



lla^ 



(5\ M+_+_ 



llg^ 
45 



-7 , 



(A.3) 



At |C|,/3,7 > 1 one has 
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M++++ = 8q;^ {-1 + (1 - 2x) \n[l/x - 1] - (x^ - x + l/2)(ln2[l/x - 1] + tt^)} , 
M+__+ = M++++(x 1/(1 - X + iO)) , M+_+_ = M++++(x ^ l/(x + iO)) , 
M++__ = M+++_ = 8q;^ . (A.4) 
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Appendix B 



In this Appendix we present the functions Xi and which appear in the result for the 
Born cross section (5.1) obtained in [6] and later rederived in [2]. 
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+ 



/5 



+ 



C-7 1 1 

2/3 /3 



1 



(/3 + 7)2 
C(/3-7)c^i 



7 C + /3 



+ 



/37 
1 



a;2 
m 

1 

\ UJi 

J m 

^ 

m 7 m 



m 



+ 



2(5 



UJ3 



\2P C7(C + /3)(/3 + 7)y m 



B{S) + 



p+1 ' (/3+7)^ 

2C7 cui 



m 



/3 + 7 m 



C(C + 7)^i I (C-/9)(C + 7) ^2 ^ 



/37 
'2C7 



+ 



m 

7-C 



C/9 
6 4 
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Here 5 — C, -\- (3 -\- ^ and the functions A, B, and 5" are defined as 



dx 



/f dx 
dx ln[l - iO + - x)] , B{u) = / —r- ^ ln[l -iQ + ux{l- x)] , 

^ ^ 

^, , 1 /■ dxvw , ( \1 — iO + V x(l ~ x)]\l — iO + w x(l — x)]) 
S [u, v,w) = - / In ' V /J V 



2 J u — vwx{l — x) [ 1 + {u + V + 'w)x{l — x) 



(B.3) 



In fact, the functions A{u) and B{u) are expressed via elementary functions and S{u, v, w) 
is expressed in terms of dilogarithmic functions [6] . 
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Appendix C 



In this Appendix we present some results obtained in [11] for the total (a) and differential 
(da/dx) cross sections of photon splitting in a screened and Coulomb potentials. The form 
factor used in the calculations was taken from [27]. 

Table C.l 

The total Born cross section a in units of Z^(nb) for a Coulomb potential (row Z = 0) and 
screened potential as a function of ui and Z. The dependence on Z is due to the atomic form 
factor. 



Z 


1 


3 


10 


30 


100 


300 


1000 


3000 





1.436 X 10- 


-3 


0.795 


12.10 


31.47 


56.5 


80.6 


107.3 


131.8 


1 


1.436 X 10" 


-3 


0.795 


12.10 


31.47 


56.4 


79.0 


97.9 


107.6 


10 


1.431 X 10" 


-3 


0.795 


12.08 


31.13 


54.1 


72.0 


84.5 


90.9 


20 


1.423 X 10' 


-3 


0.794 


12.03 


30.83 


52.98 


69.58 


81.18 


87.2 


50 


1.40 X 10" 


3 


0.790 


11.92 


30.32 


51.32 


66.32 


75.72 


81.0 


68 


1.39 X 10" 


3 


0.787 


11.85 


30.04 


50.54 


64.92 


74.22 


78.09 


92 


1.378 X 10" 


-3 


0.782 


11.78 


29.71 


49.74 


63.33 


72.07 


75.50 
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Table C.2 

The Born cross section da/dx in units of Z^{h) as a function oi x = (jJ2/u)i for different wi. The 
column Z = corresponds to the a Coulomb potential. 



z 

X 





20 


50 


82 


92 


wi = m 


0.00993 

0.0508 

0.119 

0.204 

0.296 

0.381 

0.449 

0.49 


1.78 X 10"^^ 
1.36 X 10"^^ 
9.09 X 10-13 
2.42 X 10-12 
3.88 X 10-12 
4.78 X 10-12 
5.09 X 10-12 
5.08 X 10-12 


1.59 X 10-1^ 
1.27 X 10-13 

8.84 X 10-13 
2.39 X 10-12 

3.85 X 10-12 
4.74 X 10-12 
5.06 X 10-12 
5.05 X 10-12 


1.50 X 10-1^ 
1.22 X 10-13 
8.63 X 10-13 
2.34 X 10-12 
3.79 X 10-12 
4.68 X 10-12 
4.99 X 10-12 
4.98 X 10-12 


1.44 X 10-1^ 
1.19 X 10"i3 
8.47 X 10-13 
2.31 X 10-12 
3.74 X 10-12 
4.63 X 10"i2 
4.93 X 10-12 
4.93 X 10-12 


1.42 X 10-1^ 
1.17 X 10-13 
8.41 X 10-13 
2.30 X 10-12 
3.73 X 10-12 
4.61 X 10-12 
4.92 X 10-12 
4.91 X 10-12 


uji = 10m 


0.00993 

0.0508 

0.119 

0.204 

0.296 

0.381 

0.449 

0.49 


1.84 X 10"^ 
2.51 X 10"^ 
2.76 X 10-s 
2.64 X 10-^ 
2.51 X 10-^ 
2.42 X 10"^ 
2.39 X 10"*^ 
2.38 X 10-^ 


1.79 X lO-'^ 
2.49 X 10-^ 
2.75 X 10-s 
2.63 X 10-^ 
2.49 X 10-^ 
2.41 X 10-^ 
2.37 X 10-^ 
2.37 X 10-^ 


1.76 X 10-y 

2.46 X 10"8 
2.72 X 10-8 
2.60 X 10-8 

2.47 X 10-8 
2.39 X 10"8 
2.35 X 10-8 
2.34 X 10-8 


1.73 X 10"^ 

2.44 X 10-8 
2.70 X 10-8 
2.58 X 10-8 

2.45 X 10-8 
2.37 X 10"8 
2.33 X 10-8 
2.32 X 10-8 


1.72 X 10-y 

2.43 X 10-8 
2.69 X 10-8 
2.57 X 10-8 

2.44 X 10-8 
2.36 X 10-8 
2.32 X 10-8 
2.32 X 10-8 


LOi = 100m 


0.00993 

0.0508 

0.119 

0.204 

0.296 

0.381 

0.449 

0.49 


1.65 X lO-'^ 
1.49 X 10"^ 
1.28 X lO"'^ 
1.11 X 10-'^ 
1.01 X 10-^ 
9.48 X 10-^ 
9.26 X 10-^ 
9.21 X 10-^ 


1.56 X lO-'^ 
1.40 X lO-'^ 
1.20 X lO-'^ 
1.04 X 10-^ 
9.39 X 10-8 
8.85 X 10-8 
8.64 X 10-8 
8.59 X 10-8 


1.52 X 10-'^ 
1.36 X 10"'^ 
1.16 X lO-'^ 
1.01 X 10-'^ 
9.09 X 10-8 
8.56 X 10-8 
8.36 X 10-8 
8.31 X 10-8 


1.48 X 10-'^ 
1.33 X lO"'^ 
1.13 X 10"^ 
9.83 X 10-8 
8.86 X 10"8 
8.35 X 10"8 
8.15 X 10-8 
8.10 X 10-8 


1.47 X 10-^ 
1.32 X 10-^ 
1.13 X 10""^ 
9.76 X 10-8 
8.80 X 10-8 
8.29 X 10"8 
8.09 X 10-8 
8.05 X 10-8 


LOi = 1000m 


0.00993 
0.0508 
0.119 
0.204 

0.296 
0.381 

0.449 
0. 19 


3.73 X 10"'' 
3.01 X 10"^ 
2.44 X 10-^ 
2.06 X 10-^ 
1.83 X 10"^ 
1.71 X 10"^ 
1.67 X 10"^ 
l.GG X 10"' 


2.96 X 10-'^ 
2.33 X lO-'^ 
1.85 X 10-^ 
1.54 X 10-^ 
1.36 X lO-'^ 
1.27 X 10"^ 
1.23 X lO-'^ 
1.23 X 10"' 


2.77 X 10-'' 
2.18 X 10"'^ 
1.73 X 10"'^ 
1.44 X 10-^ 
1.27 X 10-^ 
1.18 X 10"'^ 
1.15 X 10"'^ 
1.1-1 X 10"' 


2.69 X 10"'' 
2.11 X lO"'^ 
1.67 X lO"'^ 
1.39 X 10-^ 
1.22 X lO-'^ 
1.14 X lO-'^ 
1.11 X lO-'^ 
1.10 X 10"' 


2.65 X 10-^ 

2.08 X 10""^ 
1.65 X 10-"^ 
1.37 X 10-^ 
1.21 X 10-^ 
1.12 X 10""^ 

1.09 X 10-^ 
1.08 X 10""^ 
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Appendix D 



Here we present the photon sphtting amphtudes exact in the parameter Za and vahd in 
the high-energy zero- mass hmit {oj^Qi-, <^3^3 ^ rn and 62,3 <^ 1) [14]. They read 
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In the Born approximation (hnear in Za term of (D.l)) all integrals can be taken with 
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the result 
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The contribution of the region A <^ ^ to the total cross section is logarithmically amplified 
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and can be obtained also within the WW approximation. The Coulomb corrections to the 
amplitudes in this region are small compared to (D.4) and have the form 
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(D.5) 



It follows from (D.5) that in this limiting case the Coulomb correction Mjfl_ is small 



while M^\_^ and M+'l. depend only on the direction of vector A, but not on its module 
(it becomes obvious after the substitution q ^ qA). 
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